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In non-ideal plasmas, the dielectric function has to be treated beyond the random phase approx- 
imation. Correlations and well as collisions have to be included. These corrections are known as 
qq ' (dynamical) local field corrections. With the help of the Zubarev approach to linear response the- 

, ory, a relaxation time approximation is proposed leading to an interpolation scheme between static 

. local field corrections and the Drude model in the long wave length limit. The approach generalizes 

£\j ■ the Mermin approximation for the dielectric function and allows for the inclusion of a dynamical 

. | collision frequency. Exploratory calculations for a classical two-component plasma at intermediate 

■ coupling are presented. 

Q: 

On ' 

i. introduction 

cj ■ 

Many experimental observablcs in the analysis of dense plasmas are directly linked to the (longitudinal) dielectric 
I , function e(k,oj). Examples range from the reflectivity and the absorption coefficient to the pair distribution function 
and the (dynamic) structure factor [l[ . While the dielectric function for weakly coupled plasmas can be well described 
r/3 , by the random phase approximation (RPA), it is necessary to include correlations into the dielectric function to 
address the physics of strongly coupled plasmas. Corrections beyond the RPA are traditionally described by the 
Qh| so called local field corrections. For the interacting electron gas, local field corrections have been investigated in 
^2 ■ great detail since the pioneering work of Hubbard [2|. Also, approximative schemes for two-component plasmas have 
been developed Q. For general wave vectors k and frequencies to, the derived expressions tend to be very involved 
"5a ■ and tedious to calculate, see It is the objective of this communication to propose a scheme which interpolates 
£>y between the static limit co — > and the long-wave length limit k — > 0. In the course of this task, we will generalize 
i-C , an approach due to Mermin [|| and derive an approximative expression for the response function of an electron-ion 
, Ph - plasma in terms of local field corrections for the electron gas and an electron-ion collision frequency. To be specific, 
we consider a fully ionized two-component plasma of electrons and ions with temperature T and electron density 
n e . The central quantities in our description are the partial density response functions Xcc'-, where c labels the 
J> ' species, l/e{k,uf) — 1 + J2 C c' Vcc'{k)Xcc'{k, uf). Local field corrections are introduced generalizing the random phase 
approximation via 

m . 



X 



<N ; v c s c ,(k,u;) = V cc ,(k)(l-G cc ,(k,u)) 

OO : +Y / V cd (k)a-G cd (k 1 La))x d °\k,u)Vl,(k,u J ) 

o 



where V cc i (k) is the Fourier transformed potential, f2o is a normalization volume, and Xc°^ is the response function 
for the non-interacting system. For G cc > = 0, the RPA is recovered. 

II. MERMIN ANSATZ EXTENDED BY LOCAL FIELD CORRECTIONS 

Following Mermin [f| , a relaxation time approximation that obeys particle number conservation, is given by 

(M) (k \ = ( XRPA,e(fc,^ + ^)XRPA,c(fc,0) \ m 
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FIG. 1: Dynamic collision frequency for solar core conditions: n e = 6.2 x 10 25 cm 3 , T — 1.6 x 10 7 K. Various approximations 
are considered. 

where 77 is a parameter to be determined outside of the Mermin approximation. While this expression shows the 
desired Drudc-like behaviour in the long-wavelength limit allowing to identify 77 = v as a collision frequency, it fails 
to improve the static limit beyond the RPA result. Specifically, we have lim^^o x(^i w ) = XRPA,e(&,0) irrespective 
of the value of v. We rectify this shortcoming of the Mermin approach by rederiving the approximation within the 
Zubarev approach to the non-equilibrium statistical operator. Starting from the Liouville-von Neumann equation for 
the statistical operator p , we approximate the general expression with the total Hamiltonian H tot and 7/ — > 0, 

doit) i 

+ j. [H tot (t),p(t)\ = - V (p(t) - Picl (t)) , 

by a relaxation time ansatz involving the external perturbation H cx t, the intra-species interactions, and a finite 
relaxation term rj accounting for the electron-ion interaction 

doit) i 

+ ^ [H kin + V ee + V u + H eKt (t),p{t)} = -77 (p(t) - Picl (tj) ■ (2) 

Using the Zubarev technique allows to impose conserved quantities as self-consistency conditions on the relevant 
statistical operator p Te \. Proceeding along the lines presented in Q, the density response function Xcc> is then given 
in linear response by correlation functions as 

Xccf(k,u) = -/?n ; e ,.J, . ~ci\\ — • ( 3 ) 

(., .) is the Kubo product and (., .} its Laplace transform. Replacing the Kubo products by response functions, the 
extended Mermin approximation reads 

V V/ \Xee(k,uj + iri) - {iuj/ri)Xee{k,0)J 

where Xee(k,uj) is the response function of the interacting one-component electron gas. This expressions still results in 
a Drude-like form for k — > 0, while the static limit now reproduces the static local field correction, lim^^o xie M ^ {k, oS) = 

Xee(k,0). 
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FIG. 2: Imaginary part of the dielectric function as a function of the frequency to for wave vector k — 0.3 k. Parameters: 
r = 0.5, 9 = 1. Extended Mermin approach compared to other approximations. 



III. DYNAMIC COLLISION FREQUENCY 



A systematic approximation for the collision frequency in dense plasmas can be accomplished by a perturbative 
treatment of the force- force correlation function, see , 

u[oj) = — T (J ; Jo ) 

eoW p j \ / u+ir) 

Jo is the current operator, lo p i is the plasma frequency. The collision frequency can be linked to a four-particle Green's 
function. In particular, various effects such as dynamical screening and strong collisions relevant in non-ideal plasmas 
can be accounted for by partial summation of diagram sets. The net collision frequency in this so-called Gould-DcWitt 
approach is obtained as 

v(u) = ^ LB H + v t (lo) - v Boln (u) , (5) 

where i/ lb (uj) is the contribution due to loop diagrams, v 1 \uj) is the summation of ladder diagrams, and the Born 
expression has to be subtracted to avoid double counting. The interested reader is referred to Ref. Q for details. 
Here, we give the final result for the first Born approximation with respect to a dynamical screened interaction, see 



£lonm e i 3 el J 7r J 7r (oj + irj + u) 1 + u>") (— u>' — u>") 

x [Imxee(q>u' + irj)lvaxu{-q,oj" + irj) - 
lmxei(q,u' +iv)I la Xie(-q,u" +iv)] > ( 6 ) 

An adiabatic approximation with inert ions can be obtained from this expression by taking Xm(S) w ) = Xii ( ( ?) < ^( a; ) 
and Xei(q,^) = 0. We illustrate this discussion by presenting the collision frequency for a two-component plasma 
at solar core conditions n e = 6.2 x 10 25 cm~ 3 , T = 1.6 x 10 7 K, see Fig. [TJ As an example, we just compare the 
full Lenard-Balescu treatment of Eq. ([6]) with the adiabatic result indicated by m, — » oo. Also, the Born result for 
a two-component system and for the adiabatic limit are shown. Most of the features are well known such as the 
difference between the two-component Born result and the adiabatic Lenard-Balescu expression at small frequencies 
due to a different account of screening. Similar, the jump in the adiabatic Lenard-Balescu expression at the plasma 
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frequency is known to be an artifact of allowing for a undamped plasmon mode. Note, that the full calculation of 
Eq. ^ does not show such a behaviour. Instead, its overall shape is very similar to the Born approximations. The 
static limit is in accordance with a static investigation of screening in a two-component plasma of electrons and ions 
performed earlier, see [Tpj . 




frequency G) [units of Q) p| ] 

FIG. 3: Imaginary part of the dielectric function as a function of the frequency ui for wave vector k = n. Parameters: 
r = 4, 6 = 1. Extended Mermin approach compared to other approximations. 



IV. EXPLORATORY CALCULATIONS FOR A CLASSICAL TWO-COMPONENT PLASMA 

We present exploratory calculations which serve as a proof of principle taking T = 0.5 and T = 4 with 9 = 1. We 
consider an adiabatic model of interacting electrons scattering on randomly distributed but inert ions. Xee{k,uj) is 
taken for a classical OCP where the static local field corrections are related to the static structure factor S(k) via 
Gee(fc) = 1 + k 2 /n 2 (1 — 1/S(k)), k being the inverse Debye screening length. We approximate G ee {k,uj) = G ee (k). 
In later applications, this has to be tuned to more realistic expressions. Also, the collision frequency is considered in 
Born approximation with respect to a static screened potential W e i(q) ~ V^^/eRPA^, 0) , see Q, 

e f2 2 f°° 1 
Rei/H = - 2 ° 2 ° / dqq 6 W 2 l (q)Su(q)-lme RPA (q,io) . (7) 

The frequency dependence of the collision frequency is neglected, v(lo) w ^(0), to uncover the frequency dependence 
given by the Mermin approximation. Again, in order to keep things simple, we consider a uniform distribution of 
ions, i.e. Su = 1. The RPA dielectric function is taken from 

The imaginary part for the response function in extended Mermin approximation is shown in Fig. [2] for T = 0.5, k = 
0.3 k and in Fig. [3] for T = 4, k — k. For comparison, the original Mermin expression, the OCP response function, and 
the RPA are presented as well. Figure visualizes the broadening of the plasmonic excitation due to the account of 
collisions in both, the original Mermin and the extended Mermin approximation. On the other hand, for small values 
of u>, the extended Mermin approach approaches the static local field correction, as can be seen in figure [3] A similar 
situation is found for rather large values of k and T = 0.5 as in shown in Fig. 3J Here, the ideal response is given as 
well. 



V. CONCLUSIONS 



In this communication, we have proposed an interpolation scheme for the response function of a two-component 
plasma between the long-wavelength and the static limit. To this end, we combine the account of collisions via the 
Mermin ansatz with the local field description for the interacting electron gas. Thus, we obtain the broadening of 
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the reponse function due to collisions in the long-wavelength limit as well as correlations beyond RPA in the static 
limit. Exploratory calculations have shown the expected limiting behavior and indicate a flattening of the plasmon 
dispersion relation as compared to the RPA. 

Improved calculations accounting for partial degeneracy, the dynamics of the collision frequency, and dynamic 
local fields in the electronic subsystem are work in progress and subject of a forthcoming publication. In particular, 
standard approximations for dynamic local field correlations in the electron gas can easily be incorporated. 
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